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Abstract 

In this paper we considered a quantum particle moving through delute Bose-Einstein condensate 
at zero temperature. In our formulation the impurity particle interacts with the gas of uncoupled 
Bogoliubov's excitations. We constructed the perturbation theory for the Green's function of the 
impurity particle with respect to the impurity-condensate interaction employing the coherent-state 
path integral approach. The perturbative expansion for the Green's function is resumed into the 
expansion for its poles with the help of the diagrammatic technique developed in this work. The 
dispersion relation for the impurity clothed by condensate excitations is obtained and effective 
mass is evaluated beyond the Golden rule approximation. 
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I. INTRODUCTION 



Significant number of tlieoretical works were devoted to tlie quantum tlieory of super- 
fluidity on a microscopic scale. The recent interest stems from several new experiments on 
the superfluid Helium and other Bose- Einstein condensates (BEC). Rotational motion of 
molecules has been extensively studied in the superfluid helium droplets ^,6]. A unique 
properties of the measured rotational spectra pose a large number questions, such as: what 
is the collective molecule/superfluid wavefunction that describes sharp rotational states ob- 
served in experiments; what are the properties of finite systems and how is the limit of a 
bulk superfluid is reached Q]? A similar microscopic phenomena were studied in the BEC of 



Sodium atoms in magnetic traps 



22l | . A linear motion of impurities was shown to be dissipa- 



tionless for the speeds below the condensate speed of sound. A large number of theoretical 
works addressed the molecule-He droplet system using imaginary time path integral Monte- 
While in certain cases remarkable agreement with 
experimental constants was obtained , those works are strictly limited to the calculation 
of statistical properties and thus provides no real understanding of the microscopic nature 
of the dissipationless motion. The latter can only be established by considering a real time 
dynamics. A number of theoretical works considered a motion of impurity through delute 
BEC. A macroscopic particle interacting with delute BEC has been considered [3^. In this 
case the motion of particle is equivalent to the BEC in a time dependent external potential. 
This problem was treated by solving time-dependent Gross-Pitaevskii equations [j]. A mi- 
croscopic particle interacting with the delute BEC in Bogoliubov's approximation has been 
considered by several authors using general Golden rule considerations 

Q,Q,3- These 



works were based on the result of Miller et al 



2l[ | which was obtained using time-independent 



perturbation theory. The Bogoliubov's treatment has also been successfully used for the in- 
vest igation_of the force acting on the impurity particle due to the quantum fluctuations in 
BEC 



23 



2J]. Some authors treated a particle strongly interacting with Bogoliubov's BEC 



and found a possibility of self localization 15|, ll6|, ll7|| . In summary, all of the previous works 



dealt with the lowest order of the perturbation theory, either using Golden rule approxima- 
tion or considering the interaction with the fluctuations of the Gross-Pitaevskii fleld. What 
we seek is the perturbation theory that can be systematically extended to an arbitrary order. 
In present paper we use the fleld theoretical methods to develop such treatment. 
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In our recent work 



25l | we described an impurity moving through BEC as a microscopic 



particle within time-dependent perturbation theory. We introduced the Hamihonian of a 
quantum particle moving within the interacting Bose gas. No assumption was made about 
the relative mass of an impurity compared to that of the Bose particles. After the introduc- 
tion of the general Hamiltonian, the Bogoliubov's approximation was made to convert the 
Hamiltonian to the diagonal form. Then the problem was reduced to the quantum particle 
moving through the gas of non-interacting Bogoliubov's excitations. The natural way to 
compute the real time dynamics of impurity is to develop an expansion considering particle- 
BEC interaction as a perturbation. Our previous as well several works of other authors delt 
with the lowest order term that corresponds to the Golden rule limit. The higher order 
terms in direct perturbative expansion are diverging with time. Thus direct perturbation 
theory does not work beyond the Golden rule limit. 

In this work we propose a method which allows to avoid such divergence and consider 
higher orders of perturbation theory. The time evolution is written in terms of coherent state 
path integral. Because of the linear dependence of the interaction Hamiltonian on the BEC 
degrees of freedom this functional integral can be reduced to the non-Gaussian integral over 
the particle trajectories. The non-Gaussian part describing particle-BEC interaction leads 
to the formal perturbation expansion. In order to prevent the appearance of the secular 
terms, the perturbation series is resummed as an expansion for the poles of the Green's 
function with the help of the diagrammatic technique developed in this work. In this way 
we obtain the expression for the self-energy operator which is responsible for the shift of 
the Green's function pole due to the particle-BEC interaction, i.e. the dispersion relation 
for the impurity dressed by the cloud of BEC excitations. The limit of the self energy at 
small momenta allows to compute the effective mass of the particle up to an arbitrary order 
of the perturbation theory. As an example of this technique we compute the term of the 
perturbation expansion that is next to the Golden rule result and its contribution to the 
energy spectrum and effective mass of the particle. 
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II. MODEL HAMILTONIAN 



Let us start with the Hamiltonian of interacting Bose particles in secondary quantization 
representation 

P Pl,P2,P 

where U{p) is the Fourier transform of the interaction potential 

f/(p) = / U{r)e'P'dr . (2) 



The Plank constant h is set to unity here and throughout the paper. We will concentrate 
on the case of delute gas Tq ^ n~^^^ where Tq is the range of potential on which f/(r) 
differs from zero significantly and n denotes density of gas. So the Fourier transform of the 
interaction potential U{p) can be replaced by its zero component Uq = J rf/(r) which is 
connected with the length of s-scattering in first order Born approximation as following 

TTl 

as = —Uo . (3) 

Then we will consider degenerate gas at zero temperature. In this case the Hamiltonian ([T]) 
can be reduced to the diagonal form with the help of the Bogoliubov's method ^]. 

Hb = Eo + J2<p)B;B^. (4) 
p 

Here the new bosonic operators and Bp create and annihilate the collective excitations 
in BEC with the spectrum 



p2 / p2 



which has the phonon-like behavior at low momenta, i.e. e{p ^ 0) = p^J nllo/m = pc, where 
c is the speed of sound. The original particle operators and bp are connected with the 
operators of Bogoliubov's excitations Bp and Bp by following relations ^ 

bp = apBp + l3pBtp , = apB; + (3pB.p . (6) 

where the transformation coefficients are 

'^P ~ / '1 i ' ' ~ / '1 i ' ' 

V /^p - 1 V /^p - 1 

e(p) +p'^/2m + hUq 

/^P = 77 • 

nUo 
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The ground state energy of BEC is given by 

Next, we will consider a single quantum particle with mass M and momentum q interact- 
ing with the environment of Bose gas discussed above. The whole system is then described 
by following Hamiltonian 

-v2 



q 

with the particle-environment interaction 

= f E ^;+k&p<-k«q • (10) 
P,k,q 

Here the bosonic operators ctq/ctq create/annihilate the particle in state |q) and the coupling 
constant g is determined as zero Fourier component of the system-environment interaction. 
After application of the Bogoluibov's transformation to the interaction flTUl) the Hamiltonian 
of the whole system ([9]) takes the form 

H = K + Yl ^<^^ + E <p)B;B^ + E «-p«<i^p + <+P^^B^) ' (11) 
q p q.p^o 



g / Np'^ 

7p 



y 2me{p) 

The ground state energy Eq + gn is now shifted with respect to the Eq due to the system- 
condensate interaction. 

III. COHERENT-STATE PATH INTEGRAL FORMULATION OF EVOLUTION 
OPERATOR 

The main task of present work is the evaluation of the dynamical quantities such as tran- 
sition amplitudes or Greens' function of the system described the Hamiltonian in bosonic 
creation/ annihilation operator representation. In order to proceed with such kind of calcu- 



lations one can employ the coherent state path integral technique 26|, |27|] . In this section 
we will give some basic introduction into the coherent state path integral formulation of the 
dynamics of the Bose many-particle systems. 
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First, let us define the coherent state of many-particle Bose system 



l{^p}) ~ I^PU ^P25 •••5 ^Pfe5 •••) 



as an eigenstate of annihilation operator 



(12) 



(13) 



One can write the matrix element of the evolution operator in coherent state basis as the 
following functional integral 



(14) 



X exp 



Here S denotes the action depending on the trajectories Zpi^r) and z^{t) 

S{{zAt)}, {z;{t)}) = . ^ z;{t)zp{t) + z J2 z;ir)zAr) + H{{z;{r)}, {zp{r)}) , (15) 

p p 

which must be evaluated with the boundary conditions 



(16) 



^p(O) = 4 , z*p{t) = z*p , 

where z* and z'^ correspond to the bra and ket states on the left hand side of Eq. f|T4l) . 
respectively. The integration in (fT4|) is performed over all trajectories depending on time, 
and the symbol D[...] is determined as 

dz*{T)dzp{T) 



D[{z;ir)},{zpir)}] = llll 



TT 



(17) 



P T 



Here we have to note that the trajectories 2;*(r) and Zp{t) are different functions and are 
not conjugated of each other. 

Now let us concentrate on the vacuum amplitude for free Bose gas described by the 
Hamiltonian H = J2p E{p)a^ap 



(0|exp(-zt5^£;(p)a+ap)|0)= / D[{^;(r)}, {;^p(r)}]e 
p 

where 5*0 is the free action 

So = y2 d,Tzl{T)GpZp{T) , Gp -- 
„ Jo 



'iSo 



(18) 



+ E{p) 



(19) 
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Since the vacuum state is defined as |0) = |{-2p} = 0), the functional integral (fT8|) has to be 
evaluated with zero boundary conditions, i.e. -2p(0) = Zp(t) = 0. Below we will need the so 
called generating functional for the vacuum amplitude which is defined by introducing some 
auxiliary sources jpi^r) and jp(r) into the action, i.e. 



A[{jp(r)},{jp(r)}] = D[{z;{r)},{zp{r)}] 



X exp 



JO 



(20) 



Note that the sources jp(T) and jp^r) as well as the trajectories z^{t) and Zp(r) are different 
functions, so they are not conjugated. The above functional integral (l20l) can be evaluated 
by introducing the new integration variables 

Zpir) = Zp{T) + 6zp{T) , 

^;ir) = z;{r) + 6z;{T) , (21) 
where the stationary trajectories ^p(t) and ^p(r) are determined by following equations 

GpZpir) + fpir) = 0, 

G;z;{T)+jp{r) = 0. (22) 

Substituting trajectories in the form fl2Tl) into the integral fl20l) and using Eqs. fl22l) . for the 
generation functional one gets 

A[{j;(r)},{jp(r)}] = exp rfr Zp(r)jp(r) j (23) 
X I D[{5z;(r)},{5^p(r)}]exp^*rfr^5^;(r)Gp5^p(r) 

The remaining integral over deviations 6Zp{r) and 5%(t) in Eq. fl23|) exactly corresponds to 
the vacuum amplitude for free Bose particles described by normally ordered Hamiltonian and 
hence equals unity. Finally, after substituting the solution of the equation for the stationary 
trajectories fl22l) . generating functional takes the form 

A[{j;(r)},{jp(r)}] = exp dr j^rV^ir - T')jp{T)j;{r')^ , (24) 
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where Gp(t) is the Green's function of the operator Gp 

Gp(r-r') = e{T - t') exp[-iE{p){T - t')] 



(25) 



and 0(t) denotes the Heaviside step function. 

The last point is to determine the mean value of some functional of trajectories 
Ai^pi'^)}, U*pi'r)}] as follows 



-iSo 



(26) 



In the next section we will have to evaluate the mean products of the trajectories taken in 
different moments of time which can be written as the functional derivative of the generating 
functional 



^2n 



(^jpi(si)5jp2(s2)...5jpjs„)5j;^^^(s„+i)5j;^^2(s„+2)...5j;,js2„) 

xA[{j;(r)},{jp(r)}]. 



ip=jp=o 



(27) 



For example, the one particle Green's function reads 

{Zp{s)z;,{sf)) 



Gp(s-s')W 



(28) 



ip=jp=o 

Since in the case s = s' the Green's function G'(O) = (a+(s)a(s)) = 0, the expression (125 
must be written in the form G{s — s') = Q{s — s' — 0) exp[—iE{p){s — s')]. 



IV. DIAGRAMMATIC TECHNIQUE FOR THE GREEN'S FUNCTION OF THE 
RELEVANT PARTICLE 

The purpose of this section is the construction of the perturbation theory for the Green's 
function of the impurity particle Gp{t) defined as the correlation function of the creation 
and annihilation operators, i.e. Qp(t) = (ap(t)+ap(0)). Let us start with the transition 
amplitude 

w.^f = {0\B{f\e-''''\tmB, (29) 
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which describes the transition of the impurity particle from some initial state \i) to some 
final state |/) while the BEC remains in its vacuum state \0)b, i-e. state with an absence of 
Bogoliubov's excitations. It is clear that in this case \i) = |/). Below we will consider the 
eigenstate of momentum |p) as the initial and final states of the impurity 

w;.=/=|p) = = (0|B(p|e-^^*|p)|0)B . (30) 

One can write the transition amplitude as the following correlation function 

Spit) = (0|B(0|ape-'^*a+|0)|0)B = (ap(t)a+(0)) . (31) 

Here |0) denotes the state of the impurity with no particle and the brackets (...) mean 
averaging over vacuum states of the BEC and the impurity. Thus we see that the transition 
amplitude of the form (l30ll coincides with the Green's function of the impurity particle, and 
its poles determine the excitation spectrum of the particle interacting with the surrounding 
BEC. Now we can employ the path integral formulation of the matrix element of evolution 
operator developed in Sec. [Tlll 

GAt) = I D[{aliT)},{a^ir)}] a^{t)a;{0) J D[{b*^ir)} , {b^ir)}] 



X exp 



where we have determined the action of the free impurity particle 



(32) 



q 

and the action of free BEC 



Sp = ^ dral{r)G^a^ir) , = -i— + — , (33) 



= E / drbl{r)r^b^{T) , fq = -z^ + e(g) . (34) 

Since the interaction part of the action in the functional integral fl32l) has linear dependence 
on the 6*(r), 6(r)-trajectories, the BEC degrees of freedom can be immediately integrated 
out in the same manner as with the calculation of the generating functional ( 12^ . 
After eliminating the BEC from fl32|) . for the correlation function Qp one gets 

Qpit) = J D[{aliT)}, {aq(r)}] ap(t)a;(0) exp [-tSp - tSj] , (35) 
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The above integral is the non-Gaussian functional integral over impurity particle trajectories 
only, and its non-Gaussian part Sj reads 

Sj = f dr fdr' <'+q(^)«q'(^)rq(r-r')a;„_q(r')aq.(r') , (36) 

Here rq(r) represents the Green's function of the operator Fq. Now our aim is to construct 
the perturbative expansion for the path integral of the form (l35l) in powers of its non- 
Gaussian part 

G^^\t) = (ap(t)a;(0)) - {a^{t)Sja;{Q)) + ^{a^{t)Sy^m + ... + ^—J-{a^(t)S^a;{Q)){37) 
The general expression for the n-th term of the above expansion reads 

{ap{t)Sjal{0)) = / dTidT2...dTn (iri(ir2...dT:^rq^(ri - r()rq2(r2 - r2)..Tq„(r„ 

'^q'l+qi 

in ) aq'^ (ri ) {t[ ) a^>i {t[)... 

qi,qi,...,q;, qi',q^',...,q^ 

^ <i,+qn(^n)«qi.(^nXi;-qn(^n)«qi;«)«p(0)) , 

where the averaging is performed by the integration over all trajectories with the weight 
exp{—iSp). 

First, let us consider the correlation functions of the trajectories aq(r) and aq(T) in the 
integrand of Eq. (!38|) . Using general formula (127!) together with the expression for the 
generating functional for the impurity particle in the form of Eq. fl2^ where the propagator 
Gq(r) is determined as the Green's function of the operator Gq in Eq. (!33|) . for the correlator 
of the particle trajectories one gets 

(a(t)a*(ri)a(ri)a*(r{)a(r{)a*(r2)a(r2)a*(r2)a(r2)...a*(r„)a(r„)a*(r^)a(r^)a*(0)) 

= Yl - P^)G{Pl - P2)...G{P2n-l - P2n)G{P2n) • (39) 

P{{rk},H}) 

We omitted the momentum indexes in the above expression for simplicity. The sum in 
the right hand side of Eq. (1391) is performed over all permutation of the time points 
Ti, r2, Tn, r(, r2, r^. This equation which we obtained using the method of generating 
functional is in fact equivalent to the well known Wick's theorem. Next, we can substitute 
the equation (139!) into the expression (138|) and do the following: instead of permuting the 
time points in the correlator (|39|1 we will permute the time points in the product of the 
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BEC propagators r(ri — r{)r(r2 — r2)...r(r„ — r^) in the integrand in the right hand side of 
Eq. fl38|) while the particle correlator has to be taken with the fixed times r/c, r^, i.e. 

(ap(t)57a;(0)) = I" dT,dT2...dT„. C dr[dT'^...dr'^ J] V{P^ - P2)T{P2 - P^) . . .T {P^n-l - P2n) 

xGit - ri)G(n - T[)GiT[ - T2)GiT2 - r^)...G(d - r„,)G(r„ - . (40) 

Now we are able to represent each term of perturbative expansion of the Green's function 
Q{t) graphically with the help of Feynman's diagrams. Let us represent the product of the 
particle propagators by the solid lines connecting the time points 



G{t - n)G{n - r[)G{T[ - r2)...G(r„ - r;)G(r:) = • • • • • (41) 

The number of vertices equals 2n where n is the order of perturbation. The right incoming 
and left outgoing plain lines correspond to the trajectories ap(0) and ap{t) in the expression 
for the correlator fl5U]) . respectively. Thus the zeroth order term of expansion is simply given 
by 

Q^;\t) = {ap{t)a;m = Gp(t) = (42) 

Then each pair of vertices has to be connected by the BEC propagator F by the all possible 
ways is accordance with the permutation of the time points. We will represent the propa- 
gators F by wiggly lines. Let us illustrate this technique with the example of the first order 
expansion term 

G^^\t) = g^^\t)- fdr f dT'Y,Gp{t-T)Gp^^{T-T')T^{T-T')Gp{T) 
Jo Jo q 

= (43) 

Each propagator line implies the sum over its momentum, while each vertex corresponds 
to the time point and integration over it. Besides, the vertex insures the momentum con- 
servation rule, i.e. the sum of momenta of all incoming lines equals the sum of momenta of 
all outgoing lines. 

Due to the property of the Green's function G{s) ~ G(s) in the product on the right 
hand side of Eq. (HTj) we have to set t > Ti > t[ > T2 > T2 > ... > Tn > > Besides, the 
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integration over time points in the expression for the n-th perturbative term Eq. (l38ll can 
be replaced as follows 

nt i*t i*t nt nt nt 

I cLti I dT[ / dT2 / dT2... / dTn / dr^ — > 
Jo Jo Jo Jo Jo Jo 

t rri rri /.T2 i-t^^^-^ rr^ 

dTi I dT[ / dT2 dr^... dTn / dr^ . (44) 

Jo Jo Jo Jo Jo 

We can connect every pair of vertices in the diagram of n-th order by the wiggly F-lines 

by {2n)\ different ways because the number of vertices is 2n. But since only the diagrams 

with positive direction of time in each propagator r(s — s/) ,s>s' will bring the non-zero 

contribution, the whole number of all diagrams L„ of n-th order is L„ = (2n)!/2". On the 

other hand side, we have L„ = C|"C2"'~^...C|, where CJ^ = is the number of all 

possible choices of pairs of vertices to be connected by the propagator T. Next, it is clear that 

in the sum of n-th order diagrams one can meet identical graphs which can be obtained from 

each other by permuting wiggly lines. Thus one can separate all diagrams giving different 

contribution as L„ = Dnn\, where Dn denotes the number of all topologically different 

graphs, and the factor n! is due to the permutation of every pairs of vertices connected by 

wiggly line. 

Now let us write down the terms of the second {D2 = 3) and the third {D^ = 15) orders 
in diagrammatic representation in accordance with the rules established above 

W = GS^{t)+ + 



(45) 
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it) = G^^^ it) (46) 




The last step is to separate all irreducible diagrams, i.e. diagrams which can not be split 
into two diagrams of lower order by cutting one plain inner line, in all orders of perturbation. 
Let us define the self-energy graph as the sum of all irreducible blocks of each perturbation 
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+ ... . 



(47) 



Here the time points r and r' in S(r — r') denote the left and the right vertices of each 
irreducible block on the right hand side of equation, respectively. Then on can rewrite the 
expansion for the Green's function ^p(t) with the help of the self-energy operator as follows 



+ 



+ ... + 



+ 



1 



(48) 



Due to the replacement (jH]) each term of perturbation is now represented by a convolution 
of the free propagators Gp and the self-energy operators Sp, i.e. 



/ dTiGp{t - Ti) / dT2Ep{Ti - T2) 

Jo Jo 



X / c/r3Gp(r2 - Tg) / dnUpir^- r^)... dTkUp{Tk-i - n)Gp{Tk) . (49) 
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So for the Laplace transform components the equation (HHI) reads 

g^iu) = G'p(^) + G'pHrpHG'p(^) + ... + G'p(^)X'pHG'pH...X'p(^)G'pH + ... 
= G^H^jpXpR • (50) 

Here the bar denotes Laplace transform. Thus at this point we obtained the Dyson's-type 
equation for the Green's function Eq. ( l50l) . The poles of the Green's function are determined 
by the solution of the spectral equation Gp^{uj) — Sp{uj) = 0. The contribution of the self- 
energy operator into the spectral equation can be interpreted as the influence of the cloud 
of the virtual BEG excitations surrounding the particle. The solution of this equation will 
be discussed in the next section. 



V. DISPERSION RELATION AND THE ENERGY SPECTRUM OF AN IMPU- 
RITY 

In the previous section we obtained the expansion for the self-energy operator in powers 
of coupling constant 

^p(r)=<Hr) + 42)(r) + ... , (51) 

where the n-th expansion term (r) consists of all diagrams with n wiggly lines in accor- 
dance with the expansion (H71) and is proportional to g'^'"'. In accordance with Eq. (l50i) the 
Green's function yields 

^p(^) = M-T?(\ r ( \ ' (52) 
to + iE{p) - 2Jp{uj) 

and its inverse Laplace transform is given by the Fourier-Mellin integral 

1 pioo+e 

^p(^)=lT.^/ OA^V'div. (53) 

Replacing the integration variable in the above integral as a; = — ifi, for the origin of the 
Green's function one gets 

where the integration is performed over the contour shown on Fig. ([1]), and Sp{VL) = Ijp{uj = 
—iQ). First, let us concentrate on the lowest order perturbative results that can be obtained 
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by taking into account the first expansion term with one loop in the expression for the self 
energy operator Eq. fj47|) 




crW(r - r') = -C i = - $^7^Gp„k(r - r')rk(r - r') . (55) 

k 

Evaluating the Laplace transform of Eq. (!55l) and subtituting it into Eq. (150|) . for the Green's 
function in one loop approximation one gets 

1 roo —iflt 

o^Ht) = / dn , . (56) 

2-^ i-oo E{p) - 1] + Ek n-i.(p-k^.(k)+.o - ^0 
The third term in the denominator is significant only if Q is close to E{p). Thus for the 
pole of the above integral we have 

As it is well known, the real part of the pole of the Green's function determines the energy 
spectrum Ei{p) = Re(r2) while the imaginary part defines the dissipation rate. Thus the 
second term in the right hand side of Eq. ( 1571) represents the correction to the energy of 
the impurity due to the interaction with BEG while the last term describes the dissipation 
process, i.e. the energy transfer between impurity and BEG. 

Let us calculate the energy of the particle at zero momentum p = 0. After performing 
the thermodynamic limit (A^ oo, V — oo, N/V = n), i.e. replacing the sum over 
momenta Ek integral j^j^ J dk, the zero point energy can be written as 

Ei'\p = 0)=gn- -^-^ / dk ,,,,,, . (58) 

' ' ^ 167r3m } 2me{k) e{k) + k^/2m ^ ' 

In order to prevent the divergence at large momenta in the above integral one has to renor- 

malize the coupling constant g by taking into account the second order Born approximation 

for the scattering length a 

27ra / 2a /■ , \ 
9 = + — dk] , 59 

where = (1/m + 1/M)~^ is the reduced mass. The energy reexpanded in powers of 

scattering length is now finite and given as 

(-1) , , 2'Kan ( 2amc , , A , , 

Ef\p = Q) = 1 + /o(m/M) , (60) 

/ iz) = ^'^^^-^H^ + ^f^ .61) 
^{z^ - lf{z^ + 1) ■ ' 
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At this point one can consider the interesting case if the impurity has the same mass as 
the condensate particle and the scattering length a equals the scattering length for the 



interaction between Bose particles in the condensate, i.e. c = yAixan/m. Then the energy 
correction reads {Iq{z = 0) = 8/3) 



-,(1)/ „ , ,N 47ran 32 jd^n 



E '{p = 0,m = M) = 1 + W_ = . (62) 

m \ 3 V TT / 

So the energy of the resting impurity in this case coincides with the chemical potential of the 
interacting Bose gas in Bogoliubov's approximation (the expression for fiB can be obtained 
using equation for the ground state energy of BEC ([8]) and can be found, for example, in 
Ref. fl). 

The probability for the particle to stay in its initial state |p) is given by Wp = \Qp{t)\'^ and 
it decays exponentially like exp(— At), where the transition rate A is given by the imaginary 
part of the pole f2o^\ i.e. 



Ap = 2nJ2li^KE{p)-E{p-k)-e{k)) 



k 



J dk^^6{E{p) - E{p - k) - e(k)) . (63) 

The above result corresponds to the Golden rule approximation and reflects the Landau's 
criterion for the energy dissipation in BEC since the integral in the right hand side of Eq. fl63|) 
is not zero only if the momentum of the impurity particle p is more than its critical value 
Pc = Mc. 

Now we turn to the contribution of the next order perturbation term represented by a 
couple of two-loop diagrams 



a(^)(fi) = ^ ^ ^ + 

E{p - k) + E{p - V) + e(k) + e(k') - 2fi 



E- 2 2 
«7k7k' 



k,k' 



;e(k) + E(p - k) - 1] - ^0)2(e(k') + E{p - k') - ^1 - iO) 



(e(k) + e(k') + E(p - k - k') - 1^ - iO) ■ ^^^^ 
The spectral equation defining the poles of the particle Green's function is now given by 

E{p) - fio + ^4'H^o) + lofi^l^) = . (65) 
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Since last two terms in the above equation are small, we will find the solution f^o iteratively, 
i.e. in the form of expansion in powers of coupling constant 

r]?) = E{p)+^a(^\E{p)) - ^^^l^/i\E{p)) + ^a^;\E{^^^ . (66) 

The third and the fourth terms in (1661) represents the correction of the order g'^ to the result 
( !57|) . One can see that the imaginary part of the pole coming from the denominator in the 
right hand side of Eq. flMl) and defining the correction to the transition rate fl63l) does not 
contradict with the Landau's criterion. Finally, the pole with the second order contribution 
can be written as follows 

)(2)/„N _ n(i) . 9 nm I 11 



X , e(0 + e(n-2.p(l + l')A. + lm(nf(p)). 

(^e{l) -2zpl/kc^ (^e{l') -2zpV/h}j (^e(/) + e(/') + 2^11' - 2^p(I + I')/A;c) 

where e(/) = + + z),h = mc and z = m/M. The integration in the right hand 

side of Eq. ( 1671) is performed over dimensionless vectors 1, 1'. The second order contribution 
to the zero point energy does not depend on the speed of sound and diverges logarithmically 
at large momenta. 

Next, we will consider the dissipationless motion of the impurity with the momentum p 
less than the critical momentum Mc, i.e. Im(fio(p)) = 0, and expand the function ^I^q\p) 
up to second order in p. In absence of anisotropy the term linear in p disappears, and the 
energy of the impurity can be written in the form 

eP{p) = E,{p = 0) + ^ , (68) 
where the effective mass of the impurity is given by following expression 

n -1 



M^f = M 



32 r ^ 2 T 



Here we introduced new dimensionless expansion parameter 



2 / \ 2 
an I m ^ 



(69) 



9M = — , Pc = Mc. (70) 

Pc \mrj 

Two dimensionless functions Ii and I2 represent a large algebraic expressions that contain 
expansion of the integral of Eq. (167|) . The evaluation of these functions (numerical for I2) 
is shown on Fig. ([2]). 
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VI. CONCLUSION 



In this work we developed a systematic perturbation theory for the quantum propagator 
of an impurity in the degenerate BEC. The expansion of a Green's function is resummed 
as an expansion of its poles by introducing the self energy operator with the help of the 
diagrammatic technique. We demonstrate the use of this theory by computing the first two 
orders of the correction to the free propagator. In this way we obtain the energy spectrum 
and effective mass of an impurity in BEC. This theory gives access to the calculation of 
the properties of realistic systems in which the impurity-BEC interaction is not necessarily 
weak. Given the ordinary expansion developed in this work one can obtain strong coupling 
expansions by employing the variational resummation of the ordinary perturbation series 



281]. We believe that the application of the variational perturbation theory would allow 



the nonperturbative calculation of the effective mass as well as finding the regime of self- 
localization of an impurity. 
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FIG. 2: The functions that enter the expression for the effective mass of an impurity in BEC as 
given by Eq. ([69]) 
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